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^— »^ ' The rural-urban migration phenomenon is analyzed by using an agent-based computa- 

tional model. Agents are placed on lattices which dimensions varying from d = 2 up to 
^ , d = 7. The localization of the agents in the lattice define their social neighborhood (rural 

or urban) not being related to their spatial distribution. The effect of the dimension of 
' lattice is studied by analyzing the variation of the main parameters that characterizes 

the migratory process. The dynamics displays strong effects even for around one million 
*fi °f sites, in higher dimensions (d = 6, 7). 
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1. Introduction 

The rural-urban migration is a very important phenomenon which occurs in devel- 
oping economies. In our previous work, see Ref. ^ we analyzed such a phenomenon 
using an agent-based computational model, considering the rural-urban migration 
as a discrete choice problem with social interactional The migratory decision was 
modelled taking into account the pressure caused by the different earnings between 
rural and urban sectors and the neighborhood influence. The influence caused by 
neighbors was modelled like in the Ising model in a two-dimensional lattice in the 
presence of an external field. 
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Simulations in this model show some emergent properties which are consistent 
with historical data of developing countries, namely: a transitional dynamics char- 
acterized by increasing of population fraction living at urban area and per capita 
income, followed by equalization of expected wages between rural and urban sectors 
(Harris-Todaro equilibrium condition) and urban concentration. In this paper we 
extend this analysis to study the behavior toward equilibrium and its properties in 
e?-dimensional lattices, 3 < d < 7, instead of d = 2 in Ref. 



2. The benchmark model 

In this section we present the general equilibrium structure of a two-sector less 
developed economy set in Ref. H to study the rural-urban migration process 
We also review the formalization of the social interaction modelled by an Ising 
like model to as proposed in Ref. H These results will be taken to carry out the 
generalizations of subsection 13. II 



2.1. The general equilibrium structure of a rural-urban economy 

There are two productive sectors considered: urban and rural. The urban sector is 
formed by firms which are specialized in the production of manufacturated goods 
whereas the rural sector is formed by farms which produce agricultural goods. 
The production of the manufacturing sector is given by^ 

Y m = SiNS, (1) 

where N u is the amount of workers of the urban sector. £i > and a > are 
parametric constants. 

The equilibrium wage of this sector is 

w m =&K~\ (2) 

where £2 > is a parametric constant. 

The aggregated production of the rural sector is given by 

Y a =&(N-N u )*, (3) 

where TV is the total number of workers of the economic system. £3 > and cj) > 
are parametric constants. 
The equilibrium wage is 

w a =Up(N-N u )*- 1 , (4) 

where £4 > is a parametric constant. 

The terms of trade between these sectors EEI is measured by the price p: 

where p > and 7 > are parametric constants. 
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Given the equilibrium urban unemployment rate, u, as determined in Ref. ^ 
using Eqs. and one can calculate the state of the urban sector. In a similar 
manner, the rural sector has its state calculated by using Eqs. @, Q and JSJ. 
The complete deduction of the equations of state of both sectors are done in Ref. 
n Typical values for parameters mentioned above are u = 0.10, £i = 144.75, £2 = 



112.59, £3 = 500.00, U = 150.00, a = 0.70, 4> = 0.30, p = 1.00, 7 = 1.25 and 



2.2. The sectorial migration model 

The migration process was modelled by an agent-based computational model. The 
decision of migrating or not is considered taking into account the difference of 
expected wages between the sectors, called deterministic private utility, and the 
influence that indi vidua ls suffer by the group they are included in, called determin- 
istic social utilityPEEl As mentioned before, only discrete choices are allowed, so 
each worker has it state defined by 0^ 6 { — 1, +1}, where <jj = — 1 represents a rural 
worker and Oi = +1 represents an urban worker. 

As mentioned above, during the decision process each worker takes into account 
explicit and observable incentives and the influence of their social neighborhood. 
Then, the total (private and social) utility is given by 



where K > and J > are parametric constants. 

The probability that each workers reviews his/her sectorial decision is given by 
the activity a defined in Ref. 10: This parameter guarantees that only a fraction of 
the total population review their decision becoming potential migrants. The proba- 
bility that a potential migrant becomes an actual migrant is given by a cumulative 
distribution: 



where (3 > is a parametric constant which measures the heterogeneity of agents. 

The probability that a worker i migrates or not depends on the probability 
calculated in Eq. J7J. Then, the higher the value of the total utility, Eq. ©, the 
higher the probability that the worker does not change his/her sectorial decision. 

3. Properties of the transitional dynamics 

In our previous work we analyzed the transitional dynamics of the rural-urban 
migration in a two-dimensional lattice. The determination of the macrostate of the 
system in each time step can be done by using Eqs. 1)115(1 . In this section we will 
briefly review some of these results a . 



/3 = 3.00. 




(6) 



a For further details see Ref. HI 
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The main variable that characterizes the migratory process is the fraction of 
workers allocated in the urban sector, n u , also called urban share. In Fig. one 
can see the urban share as function of time. The three curves of this Figure are 
plotted for different combination of the parameters K and J of Eq. JBJ. The set 
( J = 0, K > 0) plotted in Fig. ^means that agents review their sectorial position 
taking into account only the deterministic private utility, ignoring the neighborhood 
interaction. The second case, (J > 0, K > 0), both of the effects mentioned before 
are considered in the reviewing process. In the last case, (J > 0, K = 0), only 
neighborhood influence in considered in the reviewing process what make this case 
an unrealistic one. 
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Fig. 1 . Proportion of workers at urban sector as function of time for three different set of param- 
eters J and K. Circles: (J > 0, K > 0); Squares: ( J = 0, K > 0); Crosses: ( J > 0, K = 0) 



3.1. Transitional dynamics in d- dimensions 

The simulations for (i-dimensional lattices were carried out placing workers in lat- 
tices of dimensions d = 3, 4, 5, 6 and 7. This dimensional change in the lattice 
will modify the number of nearest neighbors that each worker will have and this 
relation will depend on the dimension of the lattice: 

n b = 2d; (8) 

where rib is the number of nearest neighbors and d is the dimension of the lattice. 
In all simulations the initial urban share is n u = 0.2, what means that 20% of the 
workers are located in the urban sector. All others parameters of the system are the 
same, if not, it will be mentioned. 

In Fig. the urban share is plotted as function of time. In this figure one can 
see that the increasing of the dimension of the lattice will accelerate the migratory 
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process provoking an overshooting in the urban share as observed in developing 
economies. On the other hand, smaller dimension takes the system to reach equilib- 
rium faster than higher dimension lattices. In Eq. JHJ one can see that the variation 
in the size of the lattice will change the amount of nearest neighbors that each 
worker is connected. Therefore, as the whole parameters are kept constant to all 
simulations, the behavior of the curves in Fig. [21 are strictly related to the number 
of nearest neighbors in d-dimensional lattices. 




Fig. 2. Urban share as function of time for different dimensions of lattice (J 



The linear dimension of each lattice is L 
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An overshooting in the urban share becomes more evident when the dimension is increased. 




Fig. 3. Urban share as function of time for different dimensions of lattice for J = 0. The linear 
dimension of each lattice is i(3<j) = 100, Lujn = 32, Lr^\ = 16, £(6d) = 10, ^(7d) = 7. 
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Fig. 4. Wage of the manufacturing sector w m for different dimensions of lattice. The linear 
dimension of each lattice is £(3<j) = 100, L(4,j) = 32, = 16, £(6d) = 10, ^(7d) = 7. 



The statement above is proved in Figure [3| where all the simulations were carried 
out as in Fig.JSJbut keeping parameter J = 0. The parameter J, in Eq. ©, set null 
means that the interaction among neighbors is not being considered in the reviewing 
process done by the agents. Then, this figure clearly shows that migratory process 
is being modified due the number of neighbors, i.e., the social interaction is playing 
more important role in the decision of migrate or not. 

Figures ^] and [5] are plotting of expected urban wage w m and the rural wage w a , 
respectively. Once again, the only parameter which is changed is the dimension of 
the lattice. In these figures is shown that the curves of w m and w a have a similar 
behavior independently of the dimension of lattice. It is important to mention that 
the curves plotted for d = 7 in Figs. 0] and [5] are the only ones in which the total 
number of workers N is different from the values used in lower dimension lattices 
due to numerical limitations. Our sizes are considerably smaller than the world 
record simulated in Ref. 1111 because we do not use multi-spin coding and w a and 
w m must be reevaluated after each step. Therefore, in Fig. E3 we plotted the rural 
wage as function of time for lattice of d = 7, varying N in order to check its influence 
in the results. 

The variable r e = (1 — u)w rn /w a measures the expected wages ratio. Observ- 
ing Figures ^ and one can see that r e » 1.0, what indicates that the expected 
urban wage and the rural wage converge to same value. This is the Harris- Todaro 
equilibrium condition.^ 

The deterministic total utility of each worker is calculated by using Eq. JJJl, 
where can see two different interactions acting in the worker's decision. The first 
part of right side of Eq. @ acts like an external time dependent field and in this case 
is related to the difference of wages between rural and urban sector. The intensity of 
this field is calibrated by the constant K . The second part of right of Eq. © is the 
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Fig. 5. Wage of rural sector w a for different dimensions of lattice. The linear dimension of each 
lattice is L (3d) = 100, L (4d) = 32, L (5d) = 16, L (6d) = 10, L {7d) = 7. 
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Fig. 6. Wage of rural sector w a for different linear length of lattice for 7d. Squares: L = 5; 
Diamonds: L = 6; Triangles: L = 7; Circles: L = 8. 

interaction of a worker among his/her neighborhood and it is related to the influence 
caused by the neighbors which intensity depends on the value of the constant J. 

Therefore, in Ref. to study the influence of the variation of these parameters 
we ran several simulations where the ratio J/K and the size of the L so that we 
could check their effects in the values of the urban share. Now, we apply the same 
procedure to lattices with dimensions d > 2. 

Figured has plotted the urban share n u as function of ratio J/K for different 
lattice size L. Each figure of this set is done for different dimension from d = 2, 
top left, to d = 7, right bottom . The horizontal axis of each figure is the linear 
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Fig. 7. Urban share n u as function of the ratio J/K and the lattice size L. Lighter areas corre- 
spond to higher urban share and darker areas to lower urban share. 



dimension of the lattice L, then the number of workers is obtained by the relation 
N = L d . 

All sets of Fig. reveal a similar distribution of the equilibrium values of the 
urban share, n u , for the six different dimensions simulated. The equilibrium urban 
share is slightly different when the spatial dimension d is changed. This result is in 
agreement with Fig. |2l where one can see that d modifies the speed of the migratory 
process but it has little influence in the value of equilibrium of urban share of the 
system, what explains the similar structures seen in all sets of Fig. 
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4. Conclusion 

In this paper we analyzed the rural-urban migration process by means of an agent- 
based computational model. We extend the analysis carried out in our previous 
work to lattices with dimensions up to seven. 

The variation of the dimension of the lattices slightly modify some of the results 
found in the study of a two-dimensional lattice. To explain this, it is important to 
remember that the distribution of workers in the lattice does not mean a spatial 
distribution but the definition of the neighborhood of each individual. Therefore, 
the bigger the lattice dimension the bigger the neighborhood that each worker will 
have which is the cause of the small difference in the results. 
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